home *** CD-ROM | disk | FTP | other *** search
/ Pascal Super Library / Pascal Super Library (CW International)(1997).bin / MATH / NRPAS13 / GAMMP.PAS < prev    next >
Pascal/Delphi Source File  |  1991-04-29  |  354b  |  16 lines

  1. FUNCTION gammp(a,x: real): real;
  2. VAR
  3.    gammcf,gln: real;
  4. BEGIN
  5.    IF ((x < 0.0) OR (a <= 0.0)) THEN BEGIN
  6.       writeln('pause in GAMMP - invalid arguments'); readln
  7.    END;
  8.    IF (x < (a+1.0)) THEN BEGIN
  9.       gser(a,x,gammcf,gln);
  10.       gammp := gammcf
  11.    END ELSE BEGIN
  12.       gcf(a,x,gammcf,gln);
  13.       gammp := 1.0-gammcf
  14.    END
  15. END;
  16.